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We review the present status of the hadronic light-by-light contribution 
to muon g — 2 and critically compare recent calculations. 

1. Introduction 

The muon anomalous magnetic moment g — 2 [a^ = (g — 2)/2] has been 
measured by the E821 experiment (Muon g-2 Collaboration) at BNL with 
an impressive accuracy of 0.72 ppm [1] yielding the present world average [1] 

a^fP = 11659 208.0(6.3) x 10" 10 (1) 

with an accuracy of 0.54 ppm. New experiments [21 [3] are being designed 
to measure a M with an accuracy of at least 0.25 ppm. 

On the theory side, a lot of work has been devoted to reduce the uncer- 
tainty of the Standard Model prediction. For a recent updated discussion 
see [3] where an extensive list of references for both theoretical predictions 
and experimental results can be found. 

Here, we critically review the present status of the hadronic light-by-light 
contribution whose uncertainty will eventually become the largest theoret- 
ical error. This contribution is depicted Fig. [TJ It consists of three photon 
legs coming from the muon line connected to the external electromagnetic 
field by hadronic processes. Its contribution can be written as 

u = i V f dApi [ d * P2 I 

161 P J (2tt) 4 7 (2tt) 4 q 2 plpl(pi-m 2 )(pl-m 2 ) 

x n^ a/3 (pi,p 2 ,P3) u(p')j a (fi 4 + + m)j p u(p) (2) 



* Presented by JP. 

' On leave of absence from CAFPE and Departamento de Fi'sica Teorica y del Cosmos, 
Universidad de Granada, Campus de Fuente Nueva, E-18002 Granada, Spain. 
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Fig. 1. The hadronic light-by-light contribution to the muon g — 2. 

where q = p\ + p2 + P3 ■ To get the amplitude M. in , one needs the full 
correlator H pual3 (pi,p2,P3 — > 0), with 

n pra/3 (pi,p 2 ,p 3 ) = i 3 ^ d 4 x /" d 4 y /" dSe^ 1 '^ 2 '^ 3 ^ x 

x (0\T \v p {0)V u {x)V a (y)V^(z)\ |0) (3) 

with V fJ, (x) = qQ^q (x) and Q = diag(2, — 1, — 1)/3 the quark charges. 
The photon leg with momentum p% — > couples to the magnetic field. 

Clearly, because we have two fully independent momenta many different 
energy scales are involved in the calculation of the hadronic light-by-light 
contribution to muon g — 2. This makes it difficult to obtain the full needed 
behaviour of the correlator ([3|) from known constraints. Therefore no full 
first principles calculation exists at present. The needed results cannot be 
directly related to measurable quantities either. A first exploratory lattice 
QCD calculation has been attempted [1]. 

Using 1/N C and chiral perturbation theory (CHPT) expansion counting, 
one can organize the different contributions [5]: 

• Goldstone boson exchange contributions are order iV c and start con- 
tributing at order p e in CHPT. 

• (Constituent) quark-loop and non-Goldstone boson exchange contri- 
butions are order N c and start contributing at order p 8 in CHPT. 



Goldstone boson loop contributions are order one in 1/N C and start 
contributing at order p in CHPT. 

Non-Goldstone boson loop contributions are order one in 1/N C and 
start to contribute at order p s in CHPT. 
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The two existing full calculations, [6] and [7J are based on this classifica- 
tion. The Goldstone boson exchange contribution (GBE) was shown to be 
numerically dominant in [61 [7J after strong cancellations between the other 
contributions. [8] showed that the leading double logarithm comes from the 
GBE and was positive. [U [9] found a global sign mistake in the GBE of the 
earlier work [61 [7J which was confirmed by their authors and by |XQ|. [TT] . 

Recently, Melnikov and Vainshtein pointed out new short-distance con- 
straints on the correlator ([3]) [12], studied and extended in [13]. The authors 
of [12] constructed a model which satisfies their main new short-distance 
constraints and in which the full hadronic light-by-light contribution is given 
by GBE and axial- vector exchange contributions. Here we explicitly com- 
pare and comment on the various contributions in the different calculations. 

2. "Old" Calculations: 1995-2001 

With "old" we refer to the period 1995-2001. These calculations were 
organized according to the large N c and CHPT countings discussed above 
[5]. Notice that the CHPT counting was used just as a classification tool. We 
want to emphasize once more that the calculations in [61 [7J showed that after 
several large cancellations in the rest of the contributions, the numerically 
dominant one is the Goldstone boson exchange. Here, we discuss mainly 
the work in [6], with some comments and results from [7J and [9]. 

2.1. Pseudo-Scalar Exchange 

The pseudo-scalar exchange was saturated in [61 [7J by the Goldstone 
boson exchange. This contribution is depicted in Fig. [2] with M = vr°, r/, rf . 

Rcf. [6], used a variety of 7r°7*7* form factors 

^(pi,P2) = N c /(67r)(a/f 7r )ie^p laP2 ^(plp 2 2 ) (4) 

fulfilling as many as possible QCD constraints. A more extensive analysis 
of this form factor was done [14J finding very similar numerical results. In 
particular, the three-point form factors ^F(p 2 ,P2) used in [6] had the correct 
QCD short-distance behavior 

HQ\Q 2 ) - A/Q 2 , T(Q 2 ,0) -> B/Q 2 , (5) 

when Q 2 is Euclidean. These form factors were in agreement with available 
data including the slope at the origin as well as treating the 7r°, r] and 

1 For this one and several other contributions [6] studied thus the observance of QCD 
short-distance constraints, contrary to the often stated claim that [12] is the first 
calculation to take such constraints into account, see e.g. |15| . 
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Fig. 2. A generic meson exchange contribution to the hadronic light-by-light part 
of the muon g — 2. 



7T°, r\ and rf Exchange Contribution 


10 10 x a M 


Bijnens, Pallante and Prades [8] 


8.5 ± 1.3 


Hayakawa and Kinoshita [7J 


8.3 ± 0.6 


Knecht and Nyffeler [S] 


8.3 ± 1.2 


Melnikov and Vainshtein [12] 


11.4±1.0 



Table 1. Results for the 7r°, n and rj' exchange contributions. 



rf mixing. All form factors converged for a cutoff scale fi ~ (2 — 4) GeV 
and produced small numerical differences when plugged into the hadronic 
light-by-light contribution. 

Somewhat different ^"(pfiJ^) form factors where used in [7J [9] but the 
results agree very well (after correcting for the global sign). For comparison, 
one can find the results of [BJ [TJ [§] in Tab. Q] after adding rj and rf exchange 
contributions to the dominant it one. 

2.2. Axial- Vector Exchange 

This contribution is depicted in Fig. [2] with M = A = a§, f% and possible 
other axial- vector resonances. For this contribution one needs the ^77* 
and ^7*7* form factors. Not much is known about these but there are 
anomalous Ward identities which relate them to the -P77* and P^y*^* form 
factors. 

This contribution was not calculated in [9]. Refs. [6] and [7] used nonet 
symmetry, which is exact in the large N c limit, for the masses of the axial- 
vector resonances. Their results are shown in Tab. [2] for comparison. 
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Axial- Vector Exchange Contributions 


10 10 x 


Bijnens, Pallante and Prades [6] 
Hayakawa and Kinoshita [7J 
Melnikov and Vainshtein [12j 


0.25 ± 0.10 
0.17 ± 0.10 
2.2±0.5 



Table 2. Results for the axial-vector exchange contributions. 



T P3 




Fig. 3. Quark-loop contribution, as modeled in ENJL. 

2.3. Scalar Exchange 

This contribution is shown in Fig. [5] with M = S = ao, fo and possible 
other scalar resonances. For this contribution one needs the 6*77* and Sj*j* 
form factors. Within the extended Nambu-Jona-Lasinio (ENJL) model 
used in [6], chiral Ward identities impose relations between the constituent 
quark loop and scalar exchanges. The needed scalar form factors are also 
constrained at low energies by CHPT. Ref. [6] used nonet symmetry for the 
masses. This contribution was not included by [7J and [12]. 

In [TU] it was found the leading logarithms of the scalar exchange are 
the same as those of the pion exchange but with opposite sign. Ref. [B] 
finds that sign for the full scalar exchange contribution, obtaining 

a M (Scalar) = -(0.7 ± 0.2) • 10~ 10 . (6) 

2-4- Other contributions at leading order in 1/N C . 

This includes any other contributions that are not exchange contribu- 
tions. At short-distance, the main one is the quark-loop. At long distances 
they are often modeled as a constituent quark-loop with form-factors in 
the couplings to photons. This corresponds to the contribution shown in 
Fig. [3l Ref. [6] split up the quark momentum integration into two pieces 



G 
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A [GeV] 


0.7 


1.0 


2.0 


4.0 


10 10 x a„ 


2.2 


2.0 


1.9 


2.0 



Table 3. Sum of the short- and long-distance quark loop contributions as a function 
of the matching scale A. 



Charged Pion and Kaon Loop Contributions 


10 10 x afM 


Bijnens, Pallante and Prades (Full VMD) [6] 
Hayakawa and Kinoshita (HGS) [7J 
Melnikov and Vainshtein [12J 


-1.9 ± 0.5 
-0.45 ± 0.85 
0±1.0 



Table 4. Results for the charged and Kaon loop contributions to the hadronic 
light- by- light contribution to muon g — 2. 

by introducing an Euclidean matching scale A. At energies below A, the 
ENJL model was used to compute the quark-loop contribution while above 
A a bare (partonic) heavy quark loop of mass A was used. The latter part 
scales as 1/A 2 and mimics the high energy behavior of QCD for a massless 
quark with an IR cut-off around A. Adding these two contributions yields 
a stable result as can be seen in Tab. 

2.5. NLO in 1/N C : Goldstone Boson Loops 

At next-to-leading order (NLO) in 1/N C , the leading contribution in the 
chiral counting to the correlator in ((2|), corresponds to charged pion and 
Kaon loops which can be depicted analogously to the quark-loop in Fig. 12.41 
but with pions and Kaons running inside the loop instead. 

In [B] the needed form-factors in the 7*P + P~ and j*j*P + P~ vertices 
were studied extensively. In particular which forms were fully compati- 
ble with chiral Ward identities were studied. Full vector meson dominance 
(VMD) is one model fulfilling the known constraints. The conclusion 
reached there was that there is a large ambiguity in the momentum de- 
pendence starting at order p 6 in CHPT. Both the full VMD model of [6] 
and the hidden gauge symmetry (HGS) model of [7j satisfy the known con- 
straints. Unfortunately, this ambiguity cannot be resolved since there is 
no data for 7*7* — > ir + ir~ . Adding the charged pion and Kaon loops, the 
results obtained in [6] and [7] are listed in Tab. 01 

In view of the model dependence of this contribution, the difference 
between [6jand [7J for this contribution needs to be added linearly to the 
final uncertainty of the hadronic light-by-light contribution to ct M . 



2 Note that neither the ENJL model nor any fixed order in CHPT was used. 
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Mass Mixing 


10 10 x o M 


No OPE and Nonet Symmetry with M=1.3 GeV 
New OPE and Nonet Symmetry with M= 1.3 GeV 
New OPE and Nonet Symmetry with M= M p 
New OPE and Ideal Mixing with Experimental Masses 


0.3 
0.7 
2.8 
2.2 ± 0.5 



Table 5. Results quoted in [T^] for the pseudo- vector exchange depending of the 
/i(1285) and /i(1420) resonances mass mixing. 

3. New Short-Distance Constraints: 2003-2004 

Melnikov and Vainshtein pointed out in |12| a new short-distance con- 
straint on the correlator ([3]). This constraint is for 

(T[V^ Pl )V a (p 2 )VP(-q = - Pl -P2)]| 7 (P3 - 0)} (7) 
and follows from the OPE for two vector currents when p\ ~ p\ > > Q 2 '■ 

T\V v (pi)V a (p2)] ~ (p,/p 2 ) [qQ 2 ^q](pi + Vl) (8) 

with p = (pi — P2)j2 ~ pi ~ — p2 and Q is the light quark electrical charge 
matrix ([5]). This constraint was afterward generalized in [13] . 

The authors of [12] saturated the full correlator by exchanges. The new 
OPE constraint is satisfied by introducing a pseudo-scalar exchange with the 
vertex on the P3 side point-like. This change strongly breaks the symmetry 
between the two ends of the exchanged particle in Fig. [21 Not all OPE 
constraints on the correlator are satisfied at the same time by this model, 
but in [12] they argued that this made only a small numerical difference. 

To the pseudo-scalar exchange they added an axial- vector exchange con- 
tribution which was found to be extremely sensitive to the mixing of the 
resonances /i(1285) and /i(1420) as can be seen in Tab. taken from the 
results of [12]. The authors in [T2] took the ideal mixing result for their 
final result for a^. 

4. Comparison 

Let us now try to compare the three calculations [61 [12]. In Tab. 
the results of the leading order in 1/N C are shown. The quark loop is of the 
same order and has to be added to get the full hadronic light-by-light while 
the model used in [12] is saturated just by exchanges. In the GBE the effect 
of the new OPE in [12] is a little larger than the quark loop contributions 
of [6]. It also increases the axial- vector exchange with nonet symmetry 
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Hadronic light-by-light at 0(N C 



10 



10 



x a. 



Nonet Symmetry + Scalar Exchange [B] 
Nonet Symmetry [B] 
Nonet Symmetry [7J 
New OPE and Nonet Symmetry [12] 
New OPE and Ideal Mixing [12] 



10.2 ± 1.9 
10.9 ± 1.9 
9.4 ± 1.6 
12.1 ± 1.0 
13.6 ± 1.5 



Table 6. Full hadronic light-by-light contribution to a M at 0(N C ). The difference 
between the two results of [5] is the contribution of the scalar exchange — (0.7 ± 
0.1) • 10~ 10 . This contribution was not included in [7] and [T2"] . 

from 0.3 xlO -10 to 0.7 xl0~ 10 One thus sees a reasonable agreement in 
the comparison of the Q(N C ) results of [6j[71[T2] when using the same mass 
mixing for the axial- vectors, namely, (10. 9± 1.9, 9.4 ± 1.6, 12.1 ±1.0). 

The final differences are due to the additional increase of 1.5x 10 -10 from 
the ideal mixing in the axial vector exchange in [12] and the scalar exchange 
of -0.7xl0~ 10 in [§]. 

Let us now see what are the different predictions at NLO in 1/N C . In 
|12j , the authors studied the chiral expansion of the charged pion loop using 
the HGS model used in [7J. This model is known not to give the correct 
QCD high energy behavior in some two-point functions, in particular it 
does not fulfill Weinberg Sum Rules, see e.g. [6j. Within this model, [T2] 
showed that there is a large cancellation between the first three terms of 
an expansion of the charged pion loop contribution in powers of (m^/Mp) 2 . 
It is not clear how one should interpret this. In [6] some studies of the 
cut-off dependence of this contribution were done and the bulk of their final 
number came from fairly low energies which should be less model dependent. 
However, it is clear that there is a large model dependence in the NLO in 
1/N C contributions. But simply taking it to be (0±1)- 10~ 10 as in [T2] is 
rather drastic and certainly has an underestimated error. 

Let us now compare the results for the full hadronic light-by-light con- 
tribution to a [j, when summing all contributions. The final result quoted in 
[6], [7J and [12] can be found in Tab. [71 The apparent agreement between 
the [6] and [7J is hiding non-negligible differences which numerically almost 
compensate. There are differences in the quark loop and charged pion and 
Kaon loops and [TJdoes not include the scalar exchange. 

Comparing the results of [6] and [12], we have seen several differences 
of order 1.5 • 10 -10 , differences which are not related to the one induced by 
the new short-distance constraint introduced in [12]. These differences are 
numerically of the same order or smaller than the uncertainty quoted in [6] 
but add up as follows. The different axial- vector mass mixing account for 
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Full Hadronic Light-by-Light 


10 10 x a M 


Bijnens, Pallante and Prades |6j 
Hayakawa and Kinoshita [7] 
Melnikov and Vainshtein [12] 


8.3 ± 3.2 
8.9 ± 1.7 
13.6 ± 2.5 



Table 7. Results for the full hadronic light- by- light contribution to a^. 

— 1.5 • ICR 10 , the absence of scalar exchange in [12] accounts for — 0.7 • ICR 10 
and the absence of the NLO in 1/N C contribution in [12] accounts for — 1.9 • 
ICR 10 . These model dependent differences add up to —4.1 • ICR 10 out of the 
final —5.3 • ICR 10 difference between the results in [6] and [12] . Clearly, the 
new OPE constraint found in [12] alone does not account for the large final 
difference reading of [12] seems to suggest. 

5. Conclusions and Prospects 

At present, the only possible conclusion is that the situation of the 
hadronic light-by- light contribution to a„ is unsatisfactory. However, one 
finds a numerical agreement within roughly one sigma when comparing the 
<D(N C ) results found in [S H] , see Tab. M A new full 0(N C ) calculation 
studying the full correlator with the large N c techniques developed in [16] [TT] 
and references therein, seems feasible and desirable. 

At NLO in l/N c , one needs to control both Goldstone and non-Goldstone 
boson loop contributions. The high model dependence of the Goldstone 
boson loop is clearly visible in the different results of [6] and [7j and discussed 
in [6] and [12] . For non-Goldstone boson loops, little is known on how to 
consistently treat them, a recent attempt in another context is [18] . 

In the meanwhile, we propose as an educated guess^l 

a M = (11 ± 4) • 1(T 10 . (9) 

for the hadronic light-by- light contribution. We believe that, that this num- 
ber summarizes our present understanding of the hadronic light-by-light 
contribution to a^. One can arrive at this number in several different ways: 
the short-distance constraint and the ideal mixing for the axial-vector ex- 
change should lead to some increase of the results of [61 [7] ; the scalar ex- 
change and the pion and kaon loops are expected to lead to some decrease 
of the result of [T2] ; one can also average the leading in 1 /N c results (three 
middle results of Tab. E]). The final error remains a guess but the error in 
([9]) is chosen to include all the known uncertainties. 

3 This educated guess agrees with the one presented also at this meeting by Eduardo 
de Rafael and by one of us, JB, at DESY Theory Workshop, September 2005. 
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